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Abstract: We derive new constraints on the coupling of heavy pseudoscalar (axion-like) 
particles to photons, based on the gamma ray flux expected from the decay of these particles 
into photons. After being produced in the supernova core, these heavy axion-like particles 
would escape and a fraction of them would decay into photons before reaching the Earth. We 
have calculated the expected flux on Earth of these photons from the supernovae SN 1987A 
and Cassiopeia A and compared our results to data from the Fermi Large Area Telescope. 
This analysis provides strong constraints on the parameter space for axion-like particles. For 
a particle mass of 100 MeV, we find that the Peccei-Quinn constant, f a , must be greater 
than about 10 15 GeV. Alternatively, for f a = 10 12 GeV, we exclude the mass region between 
approximately 100 eV and 1 GeV. 
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1. Introduction 

In this work, we analyze the gamma ray flux expected on Earth from the decay of heavy 
pseudoscalar particles produced in supernovae (SNe). The analysis of this flux and comparison 
with existing experimental observations provide new limits on these particles for the mass 
range from ~ 100 eV to ~ 1 GeV. 

Pseudoscalar particles occur in many extensions of the standard model. Arguably, the 
best known pseudoscalar is the axion [|l], |2|, |3|, |J, whose existence is a natural consequence 
of the Peccei-Quinn (PQ) mechanism for solving the strong CP problem. The axion emerges 
as a pseudo-Nambu-Goldstone boson from the axial Peccei-Quinn (PQ) symmetry, Upq(1), 
spontaneously broken at the PQ (or axion) scale, f a . This scale determines the strength of 
the axion interactions, as can be seen from the Lagrangian terms, 

Ant = -igakaipk^ipk + ^-aF^F^ + (1.1) 

In Eq. ( |1 . 1| ) , g a k = C^fn^j fa is the axion-fermion coupling, with the fermion mass, and 
9a-y = C 7 a/ '(2-7T f a ) is the axion-photon coupling, where a is the fine structure constant. The 
constants Cj are model-dependent parameters typically, but not necessarily, of order 1. 

In the original axion models ]|, |], [51 ^, [7|, , the axion mass is related to the PQ constant 

via 

These standard axion models have been studied extensively in the past 30 years (see, e.g., 
ref. H for a recent review). A pseudoscalar particle which possesses some, but not all, of 
the properties of the standard axion is known as an axion-like particle. These particles have 
recently been the subject of considerable interest (see, e.g., ref. [|10| for a recent review). 

For the standard axion, experimental and astrophysical constraints exclude values of the 
PQ constant below 10 9 GeV By Eq. (L2), the axion is then very light, with m a < 10 -2 
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eV (see refs. ||, |l2| for recent discussions). However, the existence of a heavy pseudoscalar 
particle is not excluded if relation (|L^) is relaxed. Several models for these heavy- axion-like 
particles (HALPs) have been proposed in the past. Examples include models of large extra 
dimensions [|1|, [l4|, 15, 16], where heavy HALPs emerge as Kaluza-Klein (KK) modes of the 
standard light axion, models with Plank-scale induced mass terms |l7j and models of axions 



interacting with a hidden sector [18, 19, EC]]. It is, therefore, phenomenologically interesting 



to explore the axion parameter space relaxing relation (1.2). 

Observe that in some cases, for example in the models Jl9| , pC| ], the HALP is a QCD axion, 
which effectively solves the strong CP-problem. Therefore, the term axion-like particle could 
be somewhat misleading. However, in general, these heavy pseudoscalar particles are not 
responsible for the solution of the strong CP problem, and cannot be considered QCD axions. 

Rather than focusing on the specific models, in this work we keep a phenomenological 
approach, and consider a general HALP, described by two independent parameters, the mass, 
m a , and the PQ scale, f a . We assume that this particle interacts with nucleons and photons 
as described by Eq. (1.1) with CV = 1 and C 7 = 0.7. Our normalization choice for C 7 is that 
of a DFSZ axion [[/], [8| and is used in order to enable comparison with the existing literature. 
For the sake of simplicity, we neglect a possible HALP-lepton interaction. As we discuss later, 
the existence of such an interaction would not significantly change our results. 

Presently, the strongest constraints on such HALPs are derived from the analysis of 
the effects of the photons from HALPs decay on some cosmological observations such as 
the extragalactic photon background, the cosmic microwave background, and the deuterium 
abundance [21]. This argument is based on the hypothesis that the HALPs are only allowed 
to decay into photons. The analysis applies to the mass region between 1 keV and 1 GeV, 
and f a is restricted to being greater than 10 10 GeV (see Fig. 1 in ref. pi]] ) for an axion-photon 
coupling with C 7 = 0.7. 

In this paper, we extend the bound on f a to greater values for axion masses in a similar 
range, using Fermi-LAT observations of supernovae photons. 

Weakly interacting pseudoscalars may be produced via nucleon bremsstrahlung in the 
core of supernovae. Due to their weak interactions, they would freely stream out and even- 
tually decay into photons. The decay probability depends on both the HALP's mass and 
coupling with photons. In general, a fraction of the HALPs produced would decay before 
reaching Earth, generating an observable gamma-ray flux. The signal, peaked at energies of 
the order of the SN core temperature, would be stretched on a much longer timescale than 
the corresponding neutrino signal, and so the decay photons would still be observable today. 

We have analyzed the gamma-ray signal expected from the decay of massive HALPs 
produced by SN 1987A and Cassiopeia A (Cas A). We find that, for a large region of the 
axion parameter space, the gamma ray flux from the directions of SN 1987A and Cas A 
should be easily observable by the Fermi Large Area Telescope (Fermi-LAT). Observations of 
this signal would provide information on the properties of these particles. The lack of such 
observations allows us to place limits on the HALP-photon coupling (or equivalently, f a ) as a 
function of the HALP mass. We find that for f a = 10 12 GeV, the HALP mass is constrained 
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in the range from 100 eV to 1 GeV. The most stringent constraint on the Peccei-Quinn scale 
applies to HALP masses of approximately 100 MeV and restricts f a to values greater than 
about 10 15 GeV. 

This paper is organized as follows. In the next section, we analyze the HALP production 
from the SN core, calculate the corresponding photon flux from their decay, and compare our 
results to the Fermi-LAT observations, thereby obtaining limits for HALP parameters. In 
section ^, we review and discuss our results. Finally, in the appendix we provide a quantitative 
analysis of the effect of a lepton decay channel, neglected in the rest of the paper. 



2. Heavy Axion-like Particles and Supernovae 

In this section, we will first calculate the HALP luminosity from SNe, disregarding rela- 
tion ( |l.2[) , and then the flux of photons that we could expect from HALP decay on Earth. 
Combined, these calculations allow us to derive new limits on HALP properties, based on the 
possible detection of the emitted photons from HALP decay. 

2.1 Generalization of the energy loss mechanism 

Massive pseudoscalars may be produced efficiently in a SN core, due to the high temperature 
and density. With temperatures of T ~ 30 MeV and densities of p ~ 3 x 10 14 g cm -3 , the 
most efficient axion production mechanism in a SN core is the nucleon-nucleon bremsstrahlung 
process [11| 

NN^NNa, (2.1) 
where N is a nucleon. The bremsstrahlung production of massive axions has been studied in 



ref. [22], using the one pion exchange approximation. We use the same formalism to obtain 
our HALP emission rates. Setting Cm — 1 to normalize the HALP-nuclear coupling, and a 
SN core of radius of 10 km, we obtain the production rate of massive pseudoscalar in SN, 

dN ^ 1.25xl0*> (MeVs)- 1 / T \ ^ f p \ 2 x 

duo ~ /2 2 VSOMeVy V3xl0 14 gcm-3y' P 9 \2T ) ' [ ' 

where dN is the number of axions with energy between uj and uj + doo, (3 is the axion velocity 
as a fraction of the speed of light, and /12 = / a /10 12 GeV. The function <p{x) is given by 

</,(x) = ^x 2 e-*K 1 (x), (2.3) 

where K\ is the modified Bessel function of the second kind. The normalization of <j) is such 
that Jq°° (j)dx = 1. In the relativistic case ((3 ~ 1), all the energy dependence is contained in 
4>{x) which then represents the axion spectrum. The maximum of 4>{x) occurs at x — 0.62, 
approximately corresponding to an energy of 1.2 T. At high energies, 
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and thus, for m a 3> T, the production rate is suppressed by an exponential factor. This is the 
origin of the Boltzmann suppression factor expected for the production of heavy particles in 



SNe. Additionally, in the nonrelativistic regime, the production rate of Eq. (2.2) is suppressed 
by the velocity factor, /3 3 . 

As we are interested in large values of the PQ constant, corresponding to weak coupling, 
we assume the HALPs free stream from the SN core and decay into photons only when they 
are free of the SN matter. In the free-streaming regime, the total energy loss rate via HALPs, 
or HALP luminosity, is 

fdN" 



L >=L"\^r' (2 - 5) 

with dN/dw given by Eq. ( |2~2|) . 

Now we can easily generalize the novel particle cooling limit of SN 1987A, originally 
derived for a massless axion, to a general HALP. Requiring that the novel particle luminosity 
of Eq. ( |2.5| ), evaluated at typical core condition 1 s after the collapse (T ~ 30MeV, p ~ 
3 x 10 14 g cm~ 3 ) does not exceed 3 x 10 52 erg/s, we find the mass-dependent constraint on 
the PQ constant shown in Fig. ffl. Our result is in agreement with the standard bound, 
constraining f a to be greater than 10 9 GeV for masses less than 30 MeV. As expected, the 
bound is strongly reduced for larger masses. 

2.2 Constraints from decay into photons 

In the rest of this section, we discuss the stronger phenomenological constraint on axion-like 
particles which arises from their subsequent decay into photons. Weakly coupled particles 
will freely stream out of the SN core, and, if massive enough, decay into photons before 
reaching Earth, generating an observable photon flux. This flux would be observable for a 
long time after the SN explosion, as massive particles propagate more slowly than neutrinos 
and photons produced directly in the SN, and their decay time is statistically distributed. 
Therefore, a flux of gamma photons would still be observable today. 

We compare the predicted photon flux from decaying HALPs with the observations of 
gamma-rays from the Fermi Large Area Telescope, a detector designed to detect high-energy 
gamma-rays up to 300GeV, with an angular resolution of 1° or so at energies of a few lOOMeV 
& 

We study both SN 1987A and Cas A. We are interested in the photon signal from the 
decaying-HALP observable today, after a time of some 10 9 s from the direct observation of the 
SN explosion. As the SNe distances are of the order 10 12 s, decay-produced photons coming 
from an angle larger than ~ 1° will not have had time to reach Earth. Therefore, both sources 
would appear as point sources, if observed by Fermi-LAT. 

Cas A, which seems to correspond to the Flamsteed's SN of 1680, is the closest known 
core-collapse SN p3] from which a flux of gamma rays in the range we are interested in has 



not been observed pq]. Recently observed gamma rays from Cas A [26] are too energetic 
to be caused by the mechanism that we are describing here. No gamma-ray counterpart of 
SN 1987A has yet been detected. As they are not observed at the energies we are interested 
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Figure 1: Excluded region for the PQ constant as a function of the HALP mass, based on novel 
particle cooling in SN 1987A. 



in, we compare our predicted fluxes with the point source sensitivity of the Fermi-LAT, 
< 6 x lCP 9 photons cm~ 2 s -1 for photons of energies w 7 > 100 MeV, at the latitudes of the 
SNe we are considering. As we shall see below, the analysis from Cas A will provide the most 
stringent bounds on the axion-photon coupling, because of its closer distance to Earth. 

We now compute the gamma-ray flux on Earth from the decay of massive pseudoscalar 
particles produced in SNe. In the following, we assume that the HALP couples only to photons 
and nucleons. The definitions of these coupling strengths can be found after Eq. (1.1). 

The HALP lifetime is then given by 

,17 t 2 / m « r 3 



r = r a _ m ~1.7xlO i '/i2 (^J s. (2.6) 

SN 1987A occurred at a distance of c?sn — 51.4 kpc from Earth or approximately 5.3 x 10 12 s 
in natural units (c = 1). Cas A is located c?sn — 3.4 kpc or 3.5 x 10 11 s from Earth. Comparing 
these numbers with the lifetime of Eq. (|2.6| ), we see that a sizable fraction of HALPs produced 
can decay into photons before reaching Earth if their mass is large enough. The actual number 
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of HALPs that decay into photons after a time t a from production is 



dN a . 



'77 



1 



exp 



7T 



(2.7) 



where 7 = 0; /m a is the HALP Lorentz factor and lo is the HALP energy. 

In order to calculate the expected photon flux today, we need to find the HALP decay 



time t a in Eq. (|2.7|) . In natural units, the distance traversed by the HALP before decay, f5t a , 
and that by the photons produced by HALP decay which travel to Earth, t 7 , are geometrically 
related to the SN distance, dsN, via 

{I3t a ) 2 + + 2(fit a )ty cos 9 = d§ N , (2.8) 

where 9 is the photon emission angle. This can be solved for t a in terms of known parameters, 



_ (d SN + At) (1 - xp) - ^d| N (1 - x(3) 2 + /3 2 At (2d SN + At) (x 2 - 1) 
ta ~ 1 - 2x/3 + /3 2 ' (2 ' 9) 

where x = cos 9 and At = t a + 1 7 — g?sn is the HALP delay time. This is the time between the 
signal seen on Earth from a photon directly produced in the SN and that from the photons 
produced via HALP decays. In this analysis, the delay time is that between direct observation 
of each SN and present day. For SN 1987A, At ~ 22 years or 6.9 x 10 8 s and for Cas A, 
At ~ 320 years or 1.0 x 10 10 s. 

Notice that for both SN 1987A and Cas A, At <C c?sn- For the purpose of our calculation, 
we can therefore ignore the terms of order (At/dsN) 2 or higher, which contribute less than 
1% to the time t a . With this simplification t a becomes: 



1 



At | Q ( At 

- x/3 V^sn 



(2.10) 



We now proceed to calculate the resulting photon spectrum from HALPs decay. Here, we 
will make the simplifying assumption that the two photons are emitted with the same energy, 
w 7 = oj/2. This assumption is reasonable in both the nonrelativistic and the ultrarelativistic 
limits, and numerical tests show that the effect of this approximation on the calculation of 
the photon flux is small. The resulting photon spectrum is given by 



dN 1 

dbJ^ 



where 



1 



exp 



ta_\ 

It) 



dN 
doj 



f(x)dx. 



1 1-/3 2 



2 (1 - px) 



(2.11) 



(2.12) 



is the normalized probability distribution for the photon emission angle [27]. 

The expected HALP-produced photon flux on Earth can now be calculated from Eqs. 
( |2.9D , ( 2.11 ), and ( 2.12; ) in terms of the axion mass, the delay time At, and the SN distance 
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dsN- Assuming isotropic emission, massive axion-like particles would produce a flux of high 
energy photons, 

dR L = ^_dN 1 ^ 

In order to compare our prediction with the Fermi-Lat data, we need to calculate the 
expected number of photons with energy greater that 100 MeV generated by the HALPs 
decay, 

roc / dF \ 

F WQ = dw^ , (2.14) 

AooMeV \dw 7 J 

as a function of the HALP mass, m a , and the PQ constant, f a . 

Notice that since we are only interested in photons with energy greater than 100 MeV, we 
can use the approximation ( |2.4[ ) in Eq. ( |2.14 ). Setting T = 30MeV and p = 3 x 10 14 gcm~ 3 , 



and defining u>i = w 7 /lMeV and g?87A = 51.4 kpc, we find 
3/i2 2 ( ^87a\ 2 



-^100 



cm 2 s 



( ^lA f°° ^3/2^/15 ^ f U_ e -V 7 rl f(x)dx . (215) 
V «SN / JlOO J -I L J 



Using Eq. ( 2.10 ) for the axion decay time, the angular integral can be easily estimated in 
both the relativistic and nonrelativistic regimes. 

In the relativistic regime (which applies to axions with mass m a <C 200 MeV, since 
we are interested only in photons of energy above 100 MeV) the calculation of iqoo can be 
analytically performed to the end, with the use of very reasonable approximations. In this 
regime t a /jr ~ At/jr <C 1, and so the term in the square brackets in Eq. ( |2.15| ) can be 
substituted with t a /jr. In addition, we can neg lect the /3 3 ~ 1 factor in Eq. ( ^15| ). With 



these approximations, the integral in dx becomes straightforward, and Eq. ( 2.15| ) reduces to 
where = 22 yr. 

For axions heavier than 200MeV, the lower integration limit in the photon energy in 
Eq. ( 2.15| ) should be substituted with 2m a . Since this is much larger than the SN core 



temperature, the exponential in the integral dominates, and so the integral is peaked in the 
region of low energy. In this region, the exponential term is small, exp[— t a /^r] <C 1, and so 

F 100 * ^f- (^) 2 ^V**/ 15 doj, , (m a » 200MeV) . (2.17) 

The above expression cannot be explicitly calculated because of the /3 3 factor. However, since 
the integral has to be evaluated in the region where the exponential dominates, w 7 > m a ^> 
15MeV, the expression for F\qq is going to behave as 



^100 



f{ 2 2 e~ ma/i , (m a » 200 MeV), (2.18) 
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Figure 2: Flux of photons, Fioo, from HALP decay with energy greater than 100 MeV expected from 
SN Cas A (blue curve on top) and SN 1987A (red curve on the bottom) as a function of HALP mass. 
The PQ constant is f a = 10 12 GeV in both cases. The point source sensitivity of the Fermi-LAT and 
earlier gamma-ray telescope, EGRET, are also shown. 



with C ~ 30MeV, predicting an exponential suppression of the photon flux for large axion 
mass. 



We show the results of the numerical integration of Eq. ( 2.14| ) in Fig. || and |3[ In Fig. ||, we 



compare the photon flux from HALP decay, -Fiooi with the Fermi-LAT point source sensitivity 
as a function of the HALP mass, for the case of f a = 10 12 GeV. For comparison, we also show 
the detection limit from the Energetic Gamma Ray Experiment Telescope (EGRET), an older 
gamma-ray telescope with point source sensitivity of approximately 10 _7 photons cm~ 2 s _1 , for 
photon energies of o; 7 > 100 MeV, at high latitudes p8[ . Where the predictions are above the 
Fermi-LAT point source sensitivity, we are able to provide an exclusion limit. In this case, 
the region of HALP mass between 10 keV and 1 GeV is excluded from the lack of observation 
of a photon flux from SN 1987A. The data from Cas A exclude a larger mass region, from 
approximately 100 eV to lGeV. 

Fig. H shows the complete exclusion region obtained numerically, in the m a -f a plane. The 
excluded values are those for which Fermi would be able to observe our predicted spectrum 
of gamma-rays. The constraint from Cas A is more stringent, due to both a smaller value of 
dsN and a larger value of At, with respect to SN 1987A (cfr. Eq. ( 2.1 9| ) below). 



The agreement with the analytical approximation in the relativistic regime is within a 
few percent, and the exponential suppression in the very massive regime is also evident. 

The contour line in the low mass region (m a < lOOMeV) shown in Fig. || is, according to 



- 8 - 



14.5 




Figure 3: Exclusion plot for HALP models with mass up to 1 GeV, and PQ constant f a > 10 11 GeV. 
The exclusion region shows the axion parameters which would cause a photon flux greater than 6 x 
1CP 9 cm~ 2 s _1 . 



Eq. ( 2.16| ), a straight line with equation 



log (l^v)" 13 + I log 



V d SN J \At s7A ) 



1, / m a 

H — log 

2 HlMeV 



(2.19) 



The above relation gives a quantitative prediction of the effect of the SN distance and the 
delay time on the expected photon flux. 

From the figures, we see that the bounds on f a are peaked around the value m a ~ 
100 MeV, with f a restricted to being greater than approximately 10 15 GeV. The bounds are 
strongly suppressed at lower and higher masses. This is not surprising, as the the HALP flux 
is maximized if the HALP can be easily produced, but is also heavy enough to efficiently 
decay. This differs from the standard energy loss (or novel particle cooling) arguments, which 
are particularly efficient for low mass particles. 



3. Discussion and Summary 

Here we briefly summarize and comment on the key points of our analysis and discuss possible 
future directions. The main motivation for this work is the analysis of the heavy-axion-like 
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particles (HALP) parameter space, particularly in the case when the mass and PQ constant 
are not directly related. 

We have calculated the production rate of massive HALPs in a SN core and analyzed 
the photon flux expected from their decay into photons. The flux has then been compared 
to the results of the Fermi Large Area Telescope. Our main results are shown in Figs. ^ and 
||. For a HALP mass of ~ 100 MeV, we find that f a must be greater than roughly 10 15 GeV. 
For the specific case of f a = 10 12 GeV, we find that the values of the HALP mass between 
approximately 100 eV and 1 GeV are excluded. 

The idea of constraining radiatively decaying axions from gamma-ray observations of SN 
1987 A was considered already shortly after the SN 1987 A event (see, e.g., fl29| ). The original 
analysis was applied, however, only to photons from light axions decay, which reach the earth 
within a few seconds from the observation of the explosion. 

In general, the axion bounds on f a from stellar cooling are constant for masses up to 
a few T, where T is the relevant stellar temperature, and are exponentially suppressed for 
higher masses. Our bound, instead, is maximal for HALP masses around 100 MeV (a few 
times the SN core temperature of 30 MeV), but is reduced in the high and the low mass 
region. In fact, in the mechanism we have studied, two different effects play a role: HALP 
production in the SN core, and HALP decay into photons. The production rate of heavy 
HALPs is exponentially suppressed as the energy required to produce HALPs increases, but 
the decay rate into photons increases with the third power of the mass. 

We have also investigated the SN 1987A novel cooling limit that applies when the HALP 
mass is independent of the PQ constant. The standard energy loss argument from the analysis 
of the observed SN 1987A neutrino signal becomes inefficient for axion masses greater than 
the SN core temperature, approximately 30 MeV. This is quantitatively shown in Fig. |]. This 
limit is in agreement with previous bounds derived from SN 1987A for the standard axion. 

In deriving these new bounds, we have made a number of hypotheses. First, and most 
importantly, we have calculated the axion production rate in the one-pion-exchange (OPE) 
approximation, which is not completely justifiable. As shown in |3(J, this approximation 
might induce an overestimate of the production rate. Although this conclusion cannot be 
directly applied to the results in our paper, since the analysis in does not take into 
account the effects of finite axion mass, and although our bound on f a depends only on the 
fourth-root of the production rate, it is evident that our results should be considered within 
some not easily quantifiable uncertainty, as many details of the SN physics and particles 
production in its core are yet not completely known (see also the discussion in jlljj , page 
120). 

We have also assumed that our axion-like particle interacts with nucleons and photons 
only. We have not included the possibility of axion decay into leptons, even when this would 
be kinematically possible (m a > 2m e). Even if one assumes that the pseudoscalar-lepton 
coupling is absent in the Lagrangian (e.g. following hadronic axion models), the coupling 
with photons would necessarily induce a (suppressed) coupling with electrons, and so massive 
axions could have another decay channel. This effect would reduce the expected photon 
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flux and somewhat relax our bounds. Our numerical tests discussed, in the appendix, show, 
however, that this would not be a very significant effect for large PQ constant. 

Finally, we have not considered the HALP couplings to nucleons and photons as inde- 
pendent. These couplings are still related in the same manner as standard axion models, via 
the PQ constant. This may not be the case for a more generic pseudoscalar. In principle, it is 
possible to tune the axion-photon coupling to very small values without changing the axion 
coupling with nuclei |3l]] . This would, of course, modify our bounds. 

Other possible arguments could be used to constrain the region of the parameter space 
which we have analyzed. For example, for large PQ constant and small mass, axions would 
have a lifetime longer than the age of the universe. Therefore, another possible constraint 
could be derived from the analysis of the axion thermal production, and the calculation of 
the axion relic abundance today. A full calculation of the HALP relic abundance, based 
on references [32, 33], shows that this is the case for f a < 5 x 10 12 GeV, and m a > a few 
10~ 4 MeV, providing that the HALP lifetime exceeds the age of the universe. Therefore, from 
our analysis, it results that the SN argument which we used in this work is more efficient in 
constraining the region of high f a . A full analysis of other possible constraints in that region 
would, however, be very interesting and should be considered in the future. 

Our bounds provide strong implications for the models which predict heavy-axion-like 
particles, for example, the axion models in extra dimensions. A detailed analysis of these 
models is beyond the purpose of this paper, and will be the topic of future work. 
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Appendix: Inclusion of the lepton decay channel 

In our calculations we have not included the possibility of HALP decay into leptons, even 
when this would be kinematically possible {m a > 2m^). Here we want to briefly comment on 
the impact of the axion-lepton decay on our results. 

The effect of a new decay channel would be dual: first, the axion lifetime would be 



reduced, and this would affect the exponential factor in Eq. ( 2.11 ); second, the decay products 
would not just be photons. This last corresponds to a reduction of the expected photon flux 
by the factor r/r 7 , where r 7 is the HALP lifetime into photons given by Eq. (|2.6| ), and 
r = (t" 1 + T"^ 1 ) 1 is the lifetime including the lepton decay channel. For a DFSZ axion, 

6.7 x 10 10 ff 2 



n 



(m„/MeV) v 7 ! - (2m e /m a ) 



Both effects would be negligible in the case of a very massive particle. On the other hand, 
the reduction factor t/t 7 is small for small axion mass, and could potentially cause a sizable 
reduction of the photons flux coming from relativistic pseudoscalars. However, in general, 
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whenever the condition t a /^r <C 1 is satisfied, the two effects described above cancel out, as 
can be easily understood by expanding the exponential factor in Eq. ( p. 11 ). This happens 
for large values of the PQ constant, and for relatively small axion mass, as confirmed by our 
numerical analysis in Fig. |j. Therefore, we expect the reduction induced by the lepton decay 
to be sizable only in the region of intermediate axion mass, and small PQ constant. 
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Figure 4: Percentage reduction of the flux of photons with energy above 100 MeV expected on Earth 
if the lepton decay channel is open. 

Our numerical results are shown in Fig. ||[ The figure represents the percentage of the 
flux of photons with energy above 100 MeV expected on Earth if the lepton decay channel is 



open, with respect to the flux ( 2.14 ), expected if the decay into leptons is neglected. 

In the analysis we are including the decay channel into electrons and muons, and assuming 
an axiondepton coupling as for the DFSZ axion. In general, it is possible to consider models 
where the axiondepton coupling is much smaller (hadronic axion models), reducing even more 
the effects of the lepton decay channel. 

Notice that, since t a ~ At, the reduction of the flux (at large f a ) induced by the lepton 
decay channel is more pronounced for Cas A, which has a larger delay time. 



References 

[1] R. D. Peccei and H. R. Quinn, CP Conservation in the Presence of Instantons, Phys. Rev. Lett. 
38 (1977) 1440-1443. 

[2] R. D. Peccei and H. R. Quinn, Constraints Imposed by CP Conservation in the Presence of 
Instantons, Phys. Rev. D16 (1977) 1791-1797. 

[3] S. Weinberg, A New Light Boson?, Phys. Rev. Lett. 40 (1978) 223-226. 

[4] F. Wilczek, Problem of Strong P and T Invariance in the Presence of Instantons, Phys. Rev. 
Lett. 40 (1978) 279-282. 



[5] J. E. Kim, Weak Interaction Singlet and Strong CP Invariance, Phys. Rev. Lett. 43 (1979) 103. 



- 12 - 



M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Can Confinement Ensure Natural CP 
Invariance of Strong Interactions?, Nucl. Phys. B166 (1980) 493. 

A. R. Zhitnitsky, On Possible Suppression of the Axion Hadron Interactions, (in Russian), Sov. 
J. Nucl. Phys. 31 (1980) 260. 

M. Dine, W. Fischler, and M. Srcdnicki, A Simple Solution to the Strong CP Problem with a 
Harmless Axion, Phys. Lett. B104 (1981) 199. 



J. E. Kim and G. Carosi, Axions and the Strong CP Problem, arXiv : 0807 .3125 



J. Jaeckel and A. Ringwald, The Low-Energy Frontier of Particle Physics, arXiv: 1002.0329 



G. G. Raffelt, Stars as laboratories for fundamental physics: The astrophysics of neutrinos, 
axions, and other weakly interacting particles, . Chicago, USA: Univ. Pr. (1996) 664 p. 

L. D. Duffy and K. van Bibber, Axions as Dark Matter Particles, New J. Phys. 11 (2009) 
105008, |arXiv:0904.3346| . 



N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, The hierarchy problem and new dimensions 



at a millimeter, Phys. Lett. B429 (1998) 263-272, [hep-ph/9803315| 



N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, Phenomenology, astrophysics and cosmology 
of theories with sub-millimeter dimensions and TeV scale quantum gravity, Phys. Rev. D59 



(1999) 086004, ]hep-ph/9807344( ]. 

K. R. Dienes, E. Dudas, and T. Gherghetta, Invisible axions and large-radius compactifications, 



Phys. Rev. D62 (2000) 105023, friep-ph/9912455 |. 

L. Di Leila, A. Pilaftsis, G. Raffelt, and K. Zioutas, Search for solar Kaluza-Klein axions in 



theories of low- scale quantum gravity, Phys. Rev. D62 (2000) 125011, [tiep-pri/0006327|. 

Z. Berezhiani and A. Drago, Gamma ray bursts via emission of axion-like particles, Phys. Lett. 



B473 (2000) 281-290, flhep-ph/9911333 



V. A. Rubakov, Grand unification and heavy axion, JETP Lett. 65 (1997) 621-624, 
hep-ph/9703409( |. 



Z. Berezhiani, L. Gianfagna, and M. Giannotti, Strong CP problem and mirror world: The 



Weinberg-Wilczek axion revisited, Phys. Lett. B500 (2001) 286-296, ]hep-ph/0009290 |. 

L. Gianfagna, M. Giannotti, and F. Nesti, Mirror world, super symmetric axion and gamma ray 



bursts, JEEP 10 (2004) 044, |hep-ph/0409185[ |. 

E. Masso and R. Toldra, New constraints on a light spinless particle coupled to photons, Phys. 



Rev. D55 (1997) 7967-7969, piep-ph/9702275 ]. 

M. Giannotti and F. Nesti, Nucleon nucleon bremsstrahlung emission of massive axion, Phys. 



Rev. D72 (2005) 063005, |hep-ph/0505090[ ]. 

LAT Collaboration, W. B. Atwood et. ai, The Large Area Telescope on the Fermi Gamma-ray 



Space Telescope Mission, Astrophys. J. 697 (2009) 1071-1102, |arXiv : 0902 . lOgSfl . 

G. G. Pavlov, V. E. Zavlin, B. Aschenbach, J. Truemper, and D. Sanwal, The Compact Central 



Object in Cas A: A Neutron Star with Hot Polar Caps or a Black Hole?, astro-ph/9912024 



- 13 - 



S. Hannestad and G. G. Raffelt, Stringent neutron-star limits on large extra dimensions, Phys. 
Rev. Lett. 88 (2002) 071301, |hep-ph/0110067 |. 

S. Funk, Y. Uchiyama, and f. t. F. L. C. Collaboration, Fermi-LAT discovery of GeV 
gamma-ray emission from the young supernova remnant Cassiopeia A, arXiv : 1001 . 1419} 

L. Landau and E. Lifshitz, The Classical Theory of Fields, Fourth Edition: Volume 2 (Course of 
Theoretical Physics Series), . Oxford, UK: Butterworth-Hcincmann (1980) 402 p. 

EGRET Collaboration, P. Sreekumar et. ai, EGRET observations of the extragalactic gamma 
ray emission, Astrophys. J. 494 (1998) 523-534, [|astro-ph/9709257 l . 

E. W. Kolb and M. S. Turner, Limits to the Radiative Decays of Neutrinos and Axions from 
Gamma-Ray Observations of SN 1987a, Phys. Rev. Lett. 62 (1989) 509. 

C. Hanhart, D. R. Phillips, and S. Reddy, Neutrino and axion emissivities of neutron stars from 
nucleon nucleon scattering data, Phys. Lett. B499 (2001) 9-15, [ |astro-ph/ 0003445 ] . 

D. B. Kaplan, Opening the Axion Window, Nucl. Phys. B260 (1985) 215. 

M. S. Turner, Thermal Production of Not SO Invisible Axions in the Early Universe, Phys. Rev. 
Lett. 59 (1987) 2489. 

[33] E. Masso, F. Rota, and G. Zscmbinszki, On axion thermalization in the early universe, Phys. 
Rev. D66 (2002) 023004, |hep-ph/0203221 |. 



-14- 



